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Light quark masses and CKM matrix elements from lattice QCD 
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I give a brief overview of recent results from lattice QCD calculations which are relevant for the phenomenology 
of the Standard Model. I discuss, in particular, the lattice determination of light quark masses and the calculation 
of those hadronic quantities, such as semileptonic form factors, decay constants and B-parameters, which are of 
particular interest for the analysis of the CKM mixing matrix and the origin of CP violation. 



1. LIGHT QUARK MASSES 

Calculations of light quark masses are becom- 
ing, at present, one of the most important subject 
of investigation for lattice QCD. Indeed, in spite 
of their relevant phenomenological interest, light 
quark masses are still among the less known fun- 
damental parameters of the Standard Model. 

In the study of quark masses, the lattice 
method provides a unique approach: quark 
masses, defined as effective couplings renormal- 
ized at short distances, can be in fact determined 
on the lattice from non-perturbative calculations 
of hadronic quantities. An inspection of recent 
lattice results for quark masses suggests, however, 
that in such calculations a better understanding 
and quantification of the systematic errors is still 
an important requirement. To be specific, let us 
consider the lattice determinations of the strange 
quark mass, m s . A compilation of results for this 
mass, in the MS scheme, at the renormalization 
scale fi = 2 GeV, is shown in Table [j]. All the 
results have been obtained in the quenched ap- 
proximation. As can be seen from the table, al- 
though the lattice predictions all lie in the range 
100 fn s 160 MeV, however, within the errors 
quoted by the authors, the results obtained from 
different numerical studies are often in disagree- 
ment. 

The main two sources of systematic errors, 
which are responsible for the discrepancies dis- 
cussed above, are easily identified. The quark 
mass that is directly computed in lattice simu- 
lations is the (short-distance) bare lattice quark 



Table 1 

Lattice results for the strange quark mass 
m s (2GeV) in the MS scheme. We denote by 
0(a 2 ) and NPR the results obtained by using 
the non-perturbatively O (a) -improved action and 
non-perturbative renormalization respectively. 

Yr m s (2GeV)/MeV 
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128 ± 18 
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mass m(a), where a is the lattice spacing. This 
quantity is typically determined by fixing, to its 
experimental value, the mass of a hadron con- 
taining a quark with the same flavour. Since the 
non-perturbative lattice calculation of hadronic 
quantities is affected by discretization errors (in 
general of 0(a)), the same systematic errors will 
also propagate into the final determination of the 
quark masses. 

The second step in these calculations consists 
in relating the bare quark mass m(a) to the 
renormalized mass m(/x), in a given renormaliza- 
tion scheme. The connection is provided by a 
multiplicative renormalization constant, m(fi) = 
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Z m (a\i) m(a). The perturbative, renormalization 
group improved expression of Z m (a(i), at the 
NLO, has been given in ref. jy. Clearly, this 
matching introduces higher orders perturbative 
errors, which in turn affect the determination of 
the quark mass. 

It is a remarkable fact that, at present, both 
discretization and higher orders perturbative er- 
rors, in the lattice determination of light quark 
masses, can be (and in fact have been) drasti- 
cally reduced. The way to deal with discretiza- 
tion errors has been suggested by Symanzik p2| 
and consists in improving the lattice action and 
operators. This program has been realized per- 
turbatively in p^JT^ | and, more recently carried 
out at the non-perturbative level in |l5|]-[|l7j]. The 
discretization errors are reduced to 0(a^in the 
chiral limit. The residual cutoff effects are ex- 
pected to affect the determination of light quark 
masses by less than 1% pT[ . 

Higher orders perturbative errors in the evalu- 
ation of the quark mass renormalization constant 
are eliminated by adopting the non-perturbative 
renormalization prescription proposed in ref. Jl8| . 
The idea consists on imposing the renormal- 
ization conditions directly on non-perturbatively 
calculated correlation functions between external 
off-shell quark states. This prescription defines 
the mass in a specific scheme, called RI-MOM 
scheme. Lattice calculations in this scheme are 
thus completely non-perturbative. Perturbation 
theory only enters if one wants to convert the re- 
sult to the mass in the MS scheme. In this step 
(which is not necessary anyway) higher orders 
perturbative errors are expected to be negligible, 
since the relevant conversion factor is presently 
known at the N 2 LO Q. 

Among the lattice predictions for the strange 
quark mass (see Table [j]), the smallest values, 
corresponding to m s < 100 MeV, have been ob- 
tained in refs. as a result of an extrapo- 
lation to the continuum limit, a — > 0. However, 
the 0(a)-improved calculations of refs. ||,[|[ll| , 
in which discretization effects are expected to 
be negligible, predict larger values, thus casting 
some doubts on the reliability of previous extra- 
polations. 

In ref. |llj the first lattice calculation of light 



quark masses, which combined both the use of a 
non-perturbatively improved action and the non- 
perturbative renormalization technique, has been 
presented. As already noted in (To), the non- 
perturbative calculation of renormalization con- 
stants removes the discrepancies, observed in pre- 
vious calculations, between the results obtained 
by using different definitions of the bare quark 
mass (based on the vector and the axial-vector 
Ward identities respectively). It also gives values 
of the quark masses larger than those obtained by 
computing the renormalization constants in per- 
turbation theory. The relevant renormalization 
constants which enter in the lattice calculations 
of the quark masses are those of the scalar and 
pseudoscalar densities, Zs and Zp. By compar- 
ing the values of these constants obtained by us- 
ing the non-perturbative method with the predic- 
tions of one-loop (boosted) perturbation theory 
it is found that the differences can be larger than 
10% and 30% in the case of Zs and Z P respec- 
tively. Clearly, this error is introduced in the de- 
termination of the quark masses, when the renor- 
malization constants are evaluated in perturba- 
tion theory. As an example of a non-perturbative 
determination of renormalization constants, we 
show in Figure |l| the values of Zs as a function 
of the renormalization scale \i fll| . It is reassur- 
ing to find that this behaviour results in excellent 
agreement with the dependence predicted by the 
renormalization group equation. 

As a final estimate of the average light (m; = 
(m u + Wid)/2) and strange quark masses, in the 
MS scheme at the scale /i = 2 GeV, I quote the 
results of ref. fill : 



T JVLO 



™NLO 



4.9(4) MeV 
= 121(13) MeV 



and 



-N'LO 



= 4.5(4) MeV 
= 111(12) MeV 



(1) 



(2) 



at the NLO and N 2 LO respectively. In these es- 
timates discretization and higher orders pertur- 
bative errors are expected to be negligible, and 
the remaining uncertainty is mainly due to the 
quenched approximation. The values in eq. (□) 
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Figure 1. Zs, obtained by using the non- 
perturbative method, as a function of the renor- 
malization scale. The dashed and solid curves 
represent the solutions of the renormalization 
group equations at the LO and NLO respectively. 

are also in good agreement with recent results 
from QCD sum rules pc| . 

2. CKM MATRIX ELEMENTS 

In the analysis of the CKM mixing matrix, a 
particularly interesting issue is represented by the 
study of the unitarity triangle in the p-r\ plane. 
Indeed, a non trivial shape of this triangle is the 
signature of a CKM source of CP violation. 

At present, the phenomenological analysis of 
the unitarity triangle is constrained (mainly) by 
the existing measurements of B® s — B® s mixings, 
K° — K° mixing and semileptonic b — > u transi- 
tions. Once these measurements are compared 
with the theoretical predictions, constraints on 
the p-rj parameters are derived. The status of 
the art is shown in Figure from ref. [pl| . 
The theoretical predictions are based on the non- 
perturbative evaluation of some relevant hadronic 
matrix elements, which, in turn, are parameter- 
ized in terms of form factors, decay constants and 
B-parameters. The lattice method provide a re- 
liable approach to perform such calculations. In 
the following, I will present a short compilation 
of these lattice results. 

2.1. V u b from B — > ir,p semileptonic decays 

The determination of the V u b mixing angle is 
allowed by the analysis of B — > tt, p semileptonic 




Figure 2. The unitarity triangle in the p-rj 
plane and the constraints derived from Anid, 
Am s /Am d , \ ex \ and \V u b\ / \V c b\. 

decays. The reliability of lattice calculations in 
the study of heavy mesons semileptonic decays 
is provided by the analysis of D-meson decays, 
D — > K, K* Iv. In this processes the relevant mix- 
ing angle V cs is well constrained by the unitarity 
of the CKM matrix, so that the theoretical pre- 
dictions can be compared with the experimental 
measurements. The comparison is shown in Fig- 
ure ||, where the lattice results for the four rele- 
vant form factors, at zero momentum transfer, are 
presented together with their experimental values 
psft . A summary of these results is also given in 
Table 0. 



1.25 
1.00 
0.75 
0.50 
0.25 

1.25 
1.00 
0.75 
0.50 
0.25 
0.00 





' — = — * — * — 




r 


— a — 1 









V 























1.75 
1.50 



rA 1( 




=r~s: — 





—_ 1.00 
. —_ 0.75 
0.50 
-E 0.25 
H 0.00 
— 3 -0.25 



Figure 3. Lattice results for the D — > K,K*lis 
form factors, at zero momentum transfer, J2^j- 
The horizontal band indicates the present 
experimental average. 

For the B — > tt, p semileptonic decays, four dif- 
ferent lattice groups have presented results so far 
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Table 2 

Summary of lattice and experimental results for 
D — > K, K* semileptonic decays form factors at 
zero momentum transfer. The results for B — ► 



7T, p semileptonic decays, from ref. \30J , are also 
presented. 





L> -► 


K, K* 


B -> 7T,p 


Lattice 


Exp. 


Lattice 


/+(0) 


0.72(7) 


0.76(3) 


0.27(11) 


V(0) 


1.14(18) 


1.07(9) 


0.35(1°) 


M°) 


0.64(8) 


0.58(3) 


0.27(±|) 


A 2 (0) 


0.53(13) 


0.41(5) 


0.26(11) 



P^ , ^6| , p8| , ^0| . These are in quite good agreement 
one to each other. The more accurate determina- 
tion of form factors has been obtained in ref. pOj , 
from a combined analysis of _B-meson semilep- 
tonic decays and radiative B — > K*j decays. The 
results are shown in Table ||. 

2.2. Bk, the kaon bag parameter 

The bag parameter Bk is the relevant hadronic 
quantity which enters in the calculation of €k- A 
recent compilation of lattice results gives, for the 
NLO renormalization group invariant parameter, 
the estimate plj : 

B K = 0.86(6)(6) (3) 

Notice that, in present analysis of the unitarity 
triangle, four constraints are used to determine 
the values of two parameters, p and r\ (see Fig- 
ure ^). Thus, one of the constraints can also be 
removed, and the corresponding quantity can be 
determined together with p and r\. In this way, 
in ref. M] they obtain the value 

B K = 0.82±°f (4) 

in very good agreement with the lattice determi- 
nation of eq. (H) . This result also implies that sig- 
nificatively lower estimates of Bk, such as those 
obtained by using the QCD hadronic duality ap- 
proach (B K = 0.39 ± 0.10) §| or using the 
SU(3) symmetry and PCAC (B K = 1/3) @ 
are presently strongly disfavoured by the experi- 
ments. 



Table 3 

Summary of lattice results for D and B-meson 
decay constant and B parameters 



iV(MeV) 


W(MeV) 


iV(MeV) 


200(30) 


220(30) 


170(35) 


f B VB^/(MeV) 


Ib s 1 Jb 


Bb b 1 B B 


201(42) 


1.14(8) 


1.00(3) 



2.3. Amd and Am s from i?-meson decay 
constants and S-parameters 

The theoretical estimates of the -B-meson mass 
differences, Amj and Am s , can be expressed in 
terms of the pseudoscalar decay constants, fs d , , 
and _B-parameters. A recent compilation of lat- 
tice results for these quantities has beenpresented 
in ref. [Q and it is shown in Table [| Notice 
that the lattice value for fjj B , which has been 
predicted well before the first experimental mea- 
surement, it is in good agreement with the present 
experimental average, fo s = 243 ± 36 MeV |2l}| . 
The lattice results for this quantity have always 
been stable in the time. This is shown in Fig- 
ure |], where these results, obtained over a pe- 
riod of 10 years, are presented together with the 
current experimental value. The last point in 
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Figure 4. Lattice results for the pseudoscalar de- 
cay constant fn a ■ The horizontal band indicates 
the present experimental average. We refer to 
ref. J&ll for a compilation of results and corre- 
sponding references. 

the figure represents the preliminary result of 
the first C(a)-improved lattice calculation of the 
heavy mesons decay constants. It corresponds to 
f Ds = (230 ± 12) MeV ||. 

The lattice determinations of the i?-meson de- 
cay constant, fg, still represent a genuine pre- 
diction, since this constant has not been mea- 
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sured yet in the experiments. From the (overcon- 
strained) fit of the Standard Model, in ref. 
the value 



11. D. Becirevic et al, LPTHE-ORSAY-98-37. 



JbVBb = 213(21)(20) MeV 



(5) 



is obtained, in good agreement with the lattice de- 
termination (see Table |^) . An independent deter- 
mination of fs can be also derived by combining 
the lattice determination of the ratio Jb / Jd s , in 
which many of the systematic errors are expected 
to cancel, with the experimental measurement of 
fu s . Two recent lattice calculations give: 



From these values I find: 



0.77(8) 
0.75(5)(7) 



f B = (185 ±27 ±19) MeV 



APE [35] , , 

milc m [ ' 



(7) 



where the first error comes from the experimen- 
tal uncertainty on and the second one from 
the theoretical uncertainty on the ratio Jb/ Jd 3 - 
Eq. (0) is also in good agreement with the direct 
determination of fs given in Table |^. 
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